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Abstract 

Periodic frameworks with crystallographic symmetry are investigated 
from the perspective of a general deformation theory of periodic bar-and- 
joint structures in R d . It is shown that natural parametrizations provide 
affine section descriptions for families of frameworks with a specified graph 
and symmetry. A simple geometric setting for diaplacive phase transitions 
is obtained. Upper bounds are derived for the number of realizations of 
minimally rigid periodic graphs. 
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Introduction 

The notions of periodic graph and periodic framework emerged as abstractions 
of crystal structures. Inquiries about lattice sphere packings sprang from the 
same source. 

We show that a classical perspective used in the theory of positive definite 
quadratic forms and lattice sphere packins leads to natural parametrizations for 
placements of periodic graphs. In this setting, all placements with a specified 
crystallographic symmetry correspond with a certain affine section. As a re- 
sult, the deformation theory for periodic frameworks presented in our papers 
|BS2[ [BS3 is adapted in a natural way so as to encompass all cases of higher 
crystallographic symmetry. 



1 Symmetries of a periodic framework 

We adopt here definitions introduced in |BS2j . Let (G, V) be a d-periodic graph. 
The infinite graph G = (V, E) is assumed connected and when given a peri- 
odic placement (p,7r) in R d , the corresponding periodic framework is denoted 
(G, T,p, 7r). Recall that T C Aut(G) is a free Abelian group of rank d and it is 
a faithful representation of T by a lattice of translations of rank d. Moreover, 

p:V ^R d and n : V -> T{R d ) (1) 
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are related by: 

p o 7 = 77(7) o p for all 7 £ r (2) 



Relation ([2]) shows that n may be inferred from p, but most considerations about 
framework deformations and symmetries benefit from observing both functions. 
The quotient multigraph G/T is assumed to be finite and we put 

n = \V/T\ and m = \E/T\ (3) 

Periodic frameworks are abstract, idealized versions of crystalline materials and, 
like them, may posess other symmetries, besides those expressing periodicity un- 
der T. Thus, there might be a larger group of automorphisms T C S C Aut(G) 
and an extension of ir to a faithful representation of £ by a crystallographic 
group tt(£) C such that relation (J2|) would hold for all a E S. 

Considering that Aiti(G, T) is the normalizer of T in j4u£(G), a natural 
assumption for investigating this setup will be that T is normal in S, that is 
r C S C Aut(G, r). If all translational symmetries of the framework (G, T,p, ir) 
are in 7r(r), this is necessarily the case, since the subgroup of translations in a 
crystallographic group is normal. In general, the normality assumption would 
hold after replacing the initial periodicity group T by an appropriate subgroup 
of finite index r C T. Alternatively, instead of relaxing, one may refine the peri- 
odicity group by adopting all translational symmetries of the given framework. 

For these reasons, we proceed below with the study of framework symme- 
tries which correspond to graph automorphisms in the normalizer N (T) of T in 
Aut(G). Note that the quotient group N(T)/T acts naturally on the quotient 
graph G/T. It follows that N(T)/T is finite, since G/T is finite and G connected. 

Definition 1 We say that a £ N(T) = Aut(G,T) is a symmetry of the d- 
periodic framework (G, T,p, n) when the result of acting by a on the framework 
is the same as the result of acting by an isometry s G E(d), that is: 



s o p = p o a (4) 
In other words, we have a commutative diagram 



V R a 



la is (5) 



V R a 



As remarked above, it is convenient to keep p and n on an equal footing. Then 
([U becomes the equivalent, but more revealing condition: 

SOJ) = pO(T and C s O 7T = 7T O Go- (6) 
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where C s denotes the restriction to the group of translations T(R d ) of the 
conjugation by s in E(d) and C a denotes the restriction to T of the conjugation 
by a in Aut(G). 

With (p, 7r) given, it follows from ^ that s is uniquely determined by a. Indeed, 
assuming the origin in R d to be the image by p of a particular vertex vq (z G, 
we have: 

s{x) = Sx + 1, (7) 

with S = C S when R d = T{R d ) and t = s(0) = p{av a ) (8) 

Note that s is a translation if and only if C s is the identity, that is, if and only 
if a belongs to the centralizer C(T) of T in N(T): 

C(T) = {5 E N(T) \Sj = j5 for all 7 6 V} 

Thus, the set of all symmetries of (G,T,p, n) becomes a group under com- 
position. This is the symmetry group of the framework and will be denoted 
by E = E(G,r,p,7r) C N(T) C Aut(G). The periodicity group T is a nor- 
mal subgroup of E and the injective homomorphism a i->- s is an extension of 
7r : r — > T(R d ) to S — > E(d). For simplicity, this extension is also denoted 
by 7r. Since 7r(E) is a crystallographic group, we may refer to the framework 
(G, r,p, 7r) as having crystallographic symmetry E. 

2 Symmetry constraints 

For a given er £ -/V(r), we may identify all periodic placements (p, ir) for which 
a is a symmetry of the framework (G, T,p, tt). It will be convenient to give this 
description in terms of parameters based on the following choices: a complete 
set of representatives Vo, vi, v n -i for the vertex orbits of T on V and an 
isomorphism Z d w T. Later on, a complete set of edge representatives of the 
form e%j = (v, , Vj + 7,-j ) for the orbits of T on E will be implicated in obtaining 
equations for the length preservation of edges under deformation. 

Note that we allow the additive notation jv — v + 7 for the action of 7 € T on a 
vertex of the graph, which will facilitate writing the corresponding translation 
in a placement (p, w) as 7r(7)p(u) = p(v) + A, when the translation 77(7) 6 T(R d ) 
has the formula 7r(7)(a;) = a; + A and is identified with the translation vector 
A e R d . 

Besides the frequent identification T(R d ) = R d , which induces an inner product 
on the group of translations, other routine conventions and notations will 
be the following. With the chosen isomorphism Z d w T, the automorphisms 
of the free Abelian group T are given by d x d matrices with integer entries 
and determinant ±1 , that is Aut(T) is identified with the unimodular group 
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GL(d,Z). Moreover, by turning the free module Z d into the d-dimensional 
vector space 

R d = Z d ® z R^T® z R 

unimodular transformations will be conceived as linear transformations. In 
particular, conjugation by a G N(F) gives an automorphism Co- G Aut(T) and 
thereby a unimodular transformation which we denote by the same symbol 
C a G GL(d, Z). Thus, we have a representation N(F) — > GL(d, Z) with kernel 
C(T). 

When given a periodic placement (p, 7r) of (G, T), 7r gives an isomorphism 

R d ® z R^T{R d ) = R d (9) 

Thus, two bases are at play: the fixed Cartesian standard basis of R d = T{R d ) 
and the lattice basis which depends on it : Z d m T — > T(R d ). Sometimes, the 
coordinates based on the Cartesian basis are called geometric and those based 
on the periodicity lattice are called arithmetic. 

We have seen above that a necessary condition for a £ N(F) to be a symme- 
try of the framework (G,T,p,w) is that C a G GL(d,Z) should be orthogonal, 
when expressed in Cartesian coordinates. The remaining conditions refer to the 
translation part t in 0, namely 

t = p(a(v i ))-S(p(v i )) (10) 

which must be the same for alH = 0, 1, n— 1. It will be useful to express this in 
arithmetic coordinates. Let us assume that a acts on the vertex representatives 
vq,vi, v n -i according to the formulae: 

tr{v%) = ^o-(i) + 7i ( n ) 

where o~(i) is the index corresponding to the permutation effect of a on V/T and 
ji G r. Recall that the arithmetic or lattice coordines are introduced through 
7T : T -> T(R d ) = R d and the chosen identification Z d « T. Note that the 
periods Xi G F in (fTTjl will correspond with rii G and these vectors with 
integer entries depend only on a and the lattice identification Z d T. We allow 
here the same symbol C a for the conjugation given by a on T, its expression as 
an automorphism of Z d and its extension to Z d <S>z R- 

With efe, 1 < k < d the standard basis in Z d , we let = A G GL(d, R) denote 
the matrix with columns given by the lattice basis 7r(efc). We define vector 
parameters t\ by 

p{vi) = AU (12) 
and obtain from (Tl0|) and pip the following conditions: 

+ ni - C a ti = t„(j) + rij - C a tj, for 0<i,j<n-l (13) 
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Let us put u v = uj = A' A for the Gram matrix of the period lattice basis. Then, 
the orthogonality condition for C a € GL(d, Z) becomes: 

C\ toC a = LU (14) 

In summary, we have 

Proposition 2 A graph automorphism a £ N(T) = Aut(G, T) is a symmetry 
of the d-periodic framework (G, T,p, n) if and only if conditions and \1J$ 
are satisfied. 

Recall that the placement information (p, tt) enters in these equations through 
the parameters ti, < i < n — 1 and lu as described above. In subsequent 
sections, we shall elaborate on their role in describing symmetric periodic place- 
ments and symmetry preserving deformations. 

3 Parametrizations 

The deformation theory developed in |BS2j for periodic bar-and-joint frame- 
works in R d emphasized the analogy with the traditional theory of finite link- 
ages. In particular, equivalent realizations resulting from isometries applied to 
any given framework were not immediately factored out. However, enumerative 
purposes or other concerns require the quotient operation. In the finite case 
|BSlj . Cayley-Menger matrices or equivalently, Gram matrices serve the pur- 
pose. In the periodic case, crystallography and lattice theory have proven long 
ago the importance of the identification of the quotient 0(d, R)\GL(d, R) with 
the space of positive definite quadratic forms in d variables, itself represented 
by the open cone Q(d) of symmetric d x d matrices with positive eigenvalues 

[csiiGiisniw]. 

The parametrization used in the previous section follows this classical perspec- 
tive. All the information about the lattice of periods 7r(T), up to orthogonal 
transformations, is contained in the symmetric matrix lu — A* A 6 while 
the 'shift vectors' ti indicate (relative to the lattice basis) the placement of 
the vertex representatives Vi. By requesting that to — 0, equivalence under 
translation is eliminated as well. This yields 

Proposition 3 Let (G, T) be a d-periodic graph. Then all periodic placements 
in R d , up to equivalence under the group of Euclidean isometries E(d). are 
parametrized by (i? d )" -1 x £l(d), which is an open set of R dn+ ^ 2 \ 

Remarks. The vertex image sets of periodic placements are multilattices and 
this type of configurations has been considered in different contexts. While not 
implicating an edge structure, the study of multilattices envisaged in [FJ IPZj 
is related to the kinematics of phase transitions in crystalline materials. When 
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approached from the point of view of periodic sphere packings, as in [Si] |S2] . 
multilattices do acquire an edge structure from contacts between spheres. The 
resulting packing frameworks are a very particular class of periodic frameworks. 
A study of homogeneous sphere packings in dimension three has been under- 
taken by W. Fischer, E. Koch and H. Sowa in a series of papers e.g. |KSF| . The 
planar homogeneous case goes back to Niggli |N1[ IN2) . See also [F] . 

The bar-and-joint understanding of a framework brings in the (squared) 
length function for edges and the notion of deformations [BS2 . For a given d- 
periodic framework (G, ir) vertices become joints and edges become straight 
rigid bars between them. It is enough to register the (squared) length of a 
complete set of representatives for E/T. As mentioned earlier, with vertex 
representatives vq, u„_i already chosen, we may select m edge representatives 
of the form = (i)i , Vj + jij ) . For expressing the squared length in both 
geometric and arithmetic coordinates, we let € Z d stand for the vector with 
Ariij equal to the translation vector of 7r(7y). Then 

l i e tj) 2 = \\p( v 3 + Hj) -P(vi)\\ 2 = (tj + Uij - Ufojitj + n r] - U) (15) 

This gives a polynomial map with cubic (i ^ j) or linear (i = j) components. 

/ : (R d ) n - 1 x Q(d) -)• R m 

(16) 

f(ti, ...,t n -i,u) = ((tj + tiij - tifu{tj + riij - ti))ij e R m 

where the m pairs of indices ij correspond to the chosen representatives for the 
edge orbits E/T. 

The non-empty fibers of this map are configuration spaces of frameworks and 
the connected component of a framework in its configuration space gives the 
deformation space of the framework. 

For a simple comparison of this treatment with our presentation in |BS2| . 
we offer the following diagram. 

(R d ) n x GL(d, R) -£» E(d)\(R d ) n x GL(d,R) 

I II (17) 

R m 4- (i? d )"- 1 x Q(d) 

(R ) n x GL(d,R) is the parametrization used in [BS2 for periodic placements 
which include all isometric replicas of all frameworks. Hence, the fibers of the left 
vertical arrow are realization spaces for weighted periodic graphs (G,T,£), that 
is, periodic graphs with prescribed lengths for their edges. When isometries are 
factored out, we obtain, as stated above in Proposition [31 the parameter space 
(i? d ) n_1 x {1(d) with the bottom map (fTB"]) . With explicit formulae, we have the 
following description. 
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(R d ) n x GL(d, R) parametrizes periodic placements (p, n) by recording the 
positions of the n vertex representatives and the basis of the lattice of periods 
7r(r), that is, (p(vq), ...,p(v n -i),A). The left action of the isometry group E(d) 
on these parameters is given by 

u(p(v ), ...,_p(w„-i), A) = (uop(v Q ), ...,uop(D„_i), UA) (18) 

for an isometry u(x) = Ux + t, with U G 0(d, R) and t G R d . The quotient map 
q works by the formula 

q(p{v Q ), ...,p(v n -i),A) = (t 1 , ...,t„_i, u) = 

(19) 

= {A-\p{ Vl ) - p(v )), A-\p{v n ^) - p(v )), A* A) 

The left vertical arrow is the composition / o q. 

A direct enumerative consequence of the current presentation will be an 
upper bound for the number of distinct possible configurations of a minimally 
rigid periodic graph with generic edge lenght prescriptions. Recall from BS2, 
IBS3] that minimally rigid periodic graphs have m = dn+ (^j edge orbits. In the 
generic case, the corresponding edge length constraints are independent. There 
can be no more than linear constraints (i = j) among them, because all 

linear constraints affect only oj. Since the polynomial map (|16j) can be extented 
to complex projective coordinates in P m (C), we infer from Bezout's theorem the 
following bound. 

Proposition 4 Let (G, T) be a minimally rigid d-periodic graph with n — \V/T\ 
and m = \E/T\ = nd + (^j . Let fi be the number of cubic edge constraints 
7^ j) i n Then dn — d < \i < m and (G, T) has at most 3 M non-congruent 

configurations in R d for a generic prescription of edge lengths. 

Remark. This upper bound result is analogous to the one obtained in |BS1) 
for finite minimally rigid graphs. 

4 Actions and representations 

We may now return to symmetry considerations and elaborate on the affine 
nature of the symmetry constraints obtained in Proposition [2] in terms of place- 
ment parameters (t±, ...,t n -i,u>). 

Let us recall that, by definition, a G Aut{G. T) becomes a symmetry of a 
placement (p, 7f) when the effect of a on the periodic graph (G. T) is reproduced 
by the effect of an isometry s G E(d) on the image of the graph determined 
by p(V). In other words, the placements (p, 7r) and (p o er, n o C a ) must be 
equivalent under the action of E(d). Hence, in parameters (ti, t n —i, w), we 
must have one and the same point. By Proposition [21 the fixed point locus of 
such an action by a is given by an affine linear subvariety and this fact leads to 
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the obvious expectation that the action itself is expressed by an affine map in 
the parameters (ti, ...,t n -i,Lj). 

This is indeed the case, as ensuing computations will confirm. In order to 
discuss the action of Aut(G, V) on placements and the quotient parameter space 
(R d ) n ^ 1 x £!(cf) c R dn+ d) as a left action, we adopt the following convention. 

Definition 5 Let a G Aut(G, V) be an automorphism of a d-periodic graph 
(G, r). The left action of Aut(G,T) on periodic placements in R d is defined by 
the formula: 

a(p,w) = {poa- 1 ,iroC a -i) (20) 

Theorem 6 When expressed in parameters (ti, ...,t„_i,a;), the action A20\) cor- 
responds with an affine representation 

A:Aut(G,T)^Aff(dn+r\) (21) 
which factors through Aut{G, T)/T — Aut{G /T). 

For any subgroup V C £ C Aut(G,T), the periodic placements of (G,T) with 
crystallographyc symmetry £ are parametrized by the fixed locus of A(E) in 
(i? d )™- 1 x n(d), that is 

J"(£) = {xe (R d ) n - 1 x fi(d) : A{a)x = x, for all a e £} (22) 
The locus with full symmetry T(Aut(G,T)) is not empty. 

Proof: For the computation in coordinates (t\, ...,t n -i,u>), let us put 

a(ti, ...,t n -i,Lo) — (ti, ...,t n -i,&) 
and recall that to = to = 0. We have ui = A* A, with A = KC a -i , hence 

ib = {C a -,fwC a -i (23) 

Recall also that a^ 1 induces a permutation on {0, n — 1} by its effect on V/T. 
By (fTTj) and (Tl9t . we find cr~ 1 (wj) = v a - — ^0.-17^-1^) and then 

tj = C a (t a -i^ - i CT -i(o)) + {n a -i(o) - n a- 1 ( J )) (24) 

Formulae (|2"3")l and give the explicit form of the action of a, which is linear 
in the components of lj and affine in the components of tj , j = 1, ... , n — 1. 

The resulting homomorphism (f2~Tj) is obviously trivial on T. Since Aut(G, T) /T = 
Aut{G/T) is finite, the image group must have at least one fixed point (the 
barycenter of an orbit). 
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Corollary 7 There is an inclusion reversing correspondence £ i— > ^(S) be- 
tween subgroups Y C £ C Aut(G, Y) and a finite system of non-empty affine 
linear sections of (R d ) n ~ 1 x Sl(cf) which parametrize periodic placements with a 
specified crystallographic symmetry. 

It may be observed that this approach obtains periodic placements for (G, T) 
with full symmetry Aut(G, Y) realized by corresponding crystallographic groups, 
without recourse to a minimizing principle. Other methods for proving the 
existence of placements with higher symmetry rely explicitly on some 'energy 
functional' minimization technique for finding 'barycentric placements' |DF] or 
a harmonic 'standard placement' [KSj . 

5 Relaxing or refining symmetry 

Up to this point, our considerations have focused on a given d-periodic graph 
(G, r) with framework placements (G, Y,p, it) in R d . However, various problems 
may require a relaxation r C V or a refinement Y C T of the periodicity group. 
With the perspective gained in the preceding sections, we may introduce the 
following definition. 

Definition 8 Let G = (V,E) be an infinite graph and let Ti,r2 C Aut(G) 
be free Abelian groups of rank d such that the corresponding d-periodic graphs 
(G,Yi) admit periodic presentations in R d . Then, Y\ and Y2 are called com- 
mensurate when Yi D Y2 is of finite index in both groups Y\ and Y2- 

Let us observe the effect of relaxing periodicity from Y to a subgroup r C Y 
of index k. We select a complete set of representatives v Q = 0, vi, fk-x for 
T/r. Then, representatives for V/Y are given by Vi + Vj, with < i < n — 1 and 
<j < k- 1. 

When each periodicity group is identified with Z d , the inclusion of Y in Y cor- 
responds to an invertible matrix with integer entries M with det(M) = k. 
Of course, the T-periodic placements of (G, Y) are contained in the f-periodic 
placements of (G,T), and for a placement (p,ir) this inclusion takes the form: 

((ti)i,cj) >-+ ((tij)ij, Q), with t = = too (25) 

With tt(z/j) = rrij as translation vectors, we have p(vi + Vj) = p(vi) + rrij. By 
(|19P and its counterpart for Y, we find: 

A = AM, hence Q = A*A = M 1 ujM (26) 

tij = M~ 1 U + M- 1 A~ 1 m j (27) 

We already know that, from the perspective of (G, Y), the placements with 
'higher' symmetry Y are parametrized by an affine linear section and the above 
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computation confirms the expected fact that we have an affine inclusion map 
which identifies the parameter space for periodic placements of (G, T) with this 
affine linear section. It follows that the affine and convexity structure of the 
parameter space (i? d ) I1_1 x Q,(d) C R dn+ d) f or periodic placements of (G, V) is 
preserved when relaxing the lattice. 

Recall that a crystallographic group in dimension d is a discrete subgroup of 
isometries K C E(d), with a compact quotient E(d) /K. Bicberbach showed that 
the subgroup of translations in K, that is L — K n T(R d ), must be a lattice of 
rank d and is uniquely determined as the maximal free Abelian normal subgroup 
of K . Moreover, K/L is a finite group. Thus, if S C Aut(G) is isomorphic with 
a crystallographic group K C E(d), we may refer to the free Abelian normal 
subgroup r C £ corresponding to L d K, and form the d-periodic graph (G, T). 

We have seen above that, when (G, T) allows periodic placements in R d , some 
of them, namely those parametrized by .F(£), will have Euclidean symmetries 
given by some crystallographic group isomorphic with £ and K (which must 
be, according to another Bieberbach theorem, an affine conjugate of K). 

Under these circumstances, we may refer to E C Aut(G) as a crystallographic 
subgroup of Aut{G) and use the pair notation (G, £) for the graph G with 
the specified crystallographic symmetry £. It is also understood that, up to 
Euclidean isometry, the placements of (G, £) are those parametrized by .F(£). 
We note that .F(E) can be determined in the placement parameter space of any 
periodic graph (G, T) with T C T of finite index and stable under conjugation 
by £. This determination amouts to solving a linear system of equations with 
integer coefficients of the form (fT3|) and (fl4|) . corresponding to a finite set of 
transformations a 6 £ which provide generators for £/r. Verifications entirely 
similar to those performed above show that the affine and convexity structure 
of TiX>) is the same for all choices of T. 

The commensurability equivalence relation extends as follows. 

Definition 9 Two crystallographic subgroups Ei,S 2 C Aut(G) are called com- 
mensurate when Ei n £2 * s of finite index in both Si and £2- 

Remarks. It would be enough to assume Ei fl E2 of finite index in one of 
the groups, but we prefer the symmetric formulation. Clearly, in this case, the 
intersection £1 n £2 is itself a crystalographic subgroup of Aut{G). Subgroups 
of finite index in crystallographic groups can be found by simple procedures 

Relaxing or refining the symmetry of a framework are associated with certain 
variations within the framework's commensurability equivalence class. This 
language seems favorable for addressing geometric aspects of displacive phase 
transitions in crystalline materials |Bul ID] . The vertices of the infinite graph 
G may serve as labels for a subfamily or all of the atoms in some idealized 
crystal, with edges marking bonds. Under variations of temperature or pressure, 
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the same material may have phases with different crystallographic symmetry. 
Displacive phase transitions involve no bond rupture, hence the graph remains 
the same. When two phases have commensurate symmetry groups Si, £2, our 
approach gives the simplest geometrical common ground for a passage, namely 
■7-"(£i n £2), which contains both .F(£i) and -^(£2) as affine linear sections. 

Of course, as a guiding scenario, this has been formulated long ago. The new 
insight, at least at the geometric level, is that the symmetry preserving loci 
have a simple affine structure and description. However, non-linearity resurfaces 
when bonds are assumed to maintain their length. This is considered in the next 
section. 

6 Symmetry preserving deformations 

When returning to the edge squared length function (I16|) . a first simple remark 
is that / is Aut(G, Y) equivariant. To see this, it is convenient to write R m 
as the space R E / r of real-valued functions on E/Y. Then the left action of 
Aut(G, T) is simply defined as 

a(cf>) = (j) o a' 1 , for e R E/r (28) 

with the action on E and on E/Y denoted by the same symbol. Then, one easily 
verifies that 

/(<r(p,7r))=(r(/(p,7r)) (29) 

When given a framework (G, Y,p, tt) with crystallographic symmetry £ c N(Y) = 
Aut{G, r), and we want to consider only deformations which preserve this sym- 
metry (and all edge lengths), we have to restrict / to J"(£) and consider the 
fiber of f(p, tt). Since £ acts trivially on J r (£), the image by / must consist of 
points invariant under £, that is, / factors through R E ^. 

Remark. Strictly speaking, we should write T\E for the quotient by an action 
on the left. Expecting no harm, we continue with E/H and note that E/Y ^ 
E/Y, induces R E ^ -> R E / T . 

It follows that the edge length control for frameworks with crystallographic 
symmetry £ is given by a map which we allow to be denoted by the same 
symbol 

/ : T(T.) -> i? B/s (30) 

Thus, after appropriate rank computations, we obtain a setting entirely analo- 
gous to the basic case £ = T. 
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